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Abstract 

It is argued that topological disorder in amorphous solids can be described by local strains related 
to local reference crystals and local rotations. An intuitive localization criterion is formulated from 
this point of view. The Inverse Participation Ratio and the location of mobility edges in band tails 
is directly related to the character of the disorder potential in amorphous solid, the coordination 
number, the transition integral and the nodes of wave functions of the corresponding reference 
crystal. The dependence of the decay rate of band tails on temperature and static disorder are 
derived. Ab initio simulations on a-Si and experiments on a-Si:H are compared to these predictions. 

PACS numbers: 71.23.An, 71.55.Jv, 61.43. j 
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Electronic localization induced by diagonal disorder or by structural disorder has been 
intensively studied over nearly fifty years 11]. However, key properties like the energy depen- 
dence of the Inverse Participation Ratio (IPR), the location of the mobility edges and the 
decay rate of and tails are expressed in an obscure way, not directly accessible to experiment 
or simulation Perturbation theory has been applied to approximate the electron states 
of amorphous solids (AS), starting with a crystalline counterpart as zero order solution[^ 
even before Anderson's classical workj^j. In this Letter we suggest that a local formulation 
of perturbation theory is effective for the localized states confined to one distorted region 
and for the first time relate important physical quantities such as the decay rate of band 
tails and energy dependence of IPR etc. to basic material properties. 

Similar to the theory of elasticity^, the distorted regions in AS can be characterized 
by local strains referring to their local reference crystal (LRC) and local rotations. By a 
suitable choice of origin and orientation of LRC, the atomic displacements of a distorted 
region of AS relative to its LRC are small. Thus the relative change in potential energy 
for each distorted region in AS is small. Perturbation theory is justified for each distorted 
region. The semi-classical approximation (SCA)j^ can further simplify the calculation of 
scattering waves caused by a distorted region, since the de Broglie wavelength for low-lying 
excitations is of order one bond length (:^2.35 A in a-Si[6j), a distance much shorter than the 
characteristic range in . The motion of electronic 

packet under extra force of AS relative to LRC can 

We first formulate an intuitive localization criterion for the states confined to one distorted 
region. Then the IPR, the position of mobility edge and Urbach energy are related to 
the distortion relative to the LRC, the coordination number and the inter-cell transition 
integral. The predictions are consistent with available experiments. We also performed ab 
initio local density approximation (LDA) 9] and tight binding approximation (TBA) 10. 
computations on a-Si to verify our results. 

Consider a distorted region V, with linear size L. Using the primitive cell of LRC num- 
bering the atoms in V, the x-component of extra force suffered by an electron relative to 
that of LRC is 

later its characteristic value of is denoted as F. n is lattice index, Rn and are the 
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position vector and the (3th component of the static displacement of the atom n respectively. 
U{r — Rn) is the potential energy felt by an electron at r from the atom at Rn. 

A Bloch wave ip'^^^ of LRC passes through V, and in SCa[^, the change in the x compo- 
nent of the wave vector after scattering is 

Ai. ~ ^ (2) 

FL measures the magnitude of random potential in V. The phase shift 5„k of state tp^^ is 
determined by the change in momentum and the propagation path of the Bloch wave 



(3) 



IVk-Enkl 

where E^k is dispersion relation of the nth energy band of the LRC. FL"^ is the strength 
of a potential well (the product of the depth of potential well and the range of force) in 
standard scattering theory [l^. If the first coordination shell around an atom is spherically 
symmetric, the dispersion relation in TEA isjl2| 

Enk ~ EnO - Zin COS k^tt . (4) 

Here Eno is the middle of the nth band {k^a = vr/2), z is the coordination number of a cell, 
In is the transition integral for the nth band, a is the lattice constant in LRC. For a semi- 
quantitative discussion, crude dispersion relation (jlj) will not invalidate essential points. If 
the phase shift 5rik of the secondary scattering waves relative to the primary wave is ~ vr, 
then outside V, scattering waves will interfere destructively with the primary Bloch state. 
No probability amplitude appears outside V. A localized state is therefore formed inside T> 
due to the constructive interference of a Bloch state ip^^ and its secondary scattering waves. 

Bloch states of LRC at top of valence and at bottom of the conduction edges are sus- 
ceptible to the random potential. The former is shorter wave, sensitive to details of atomic 
displacements of a distorted region. The latter is long wave: a small random potential will 
easily produce a change in momentum comparable to Kk itself. In other words, states with 
small group velocity are easily localized. The group velocity of an electron in state ip'^^^ in 
TBA is vf^y. ~ sin near to bottom {k^a ~ 0) and states near to top {k^a ~ vr 

) have small f^j^. According to Eq.Q, they are more easily localized than the states in 
the middle of a band for a given random potential. For k close to ^, with TBA dispersion 
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relation (jl]), group velocity of state ipl is vl = Eq = Eq + zly is the top of the 

valence band. By Eq. ([3]), under TBA, for a valence state ipl with energy Ek, the change in 
phase shift with energy is ^ = ^^^^_^_J^^%^j^y/2 ■ For a given distorted region, Bloch states 
close to Eq will suffer larger phase shift. They are more readily localized than the states 
in the middle of the band. Similar conclusion holds for the Bloch states in the bottom of 
conduction band. In Fig. [T] the IPR is plotted against electron energy for a model of a-Si. 
Large IPR appears at the edges of a band, in agreement with the above prediction. 

The upper mobility edge of the valence band is the deepest energy level ii^^^that the 
largest distorted region could localize, i.e. produce a phase shift vr for the corresponding 
Bloch state. In TBA, this leads to sin k^a = The energy difference between the top 

of a band and the mobility edge is E^^ = zyl^ {I - [1 - {-^^jv^fY''^} ~ -^^^ last ~ only 
holds for ^^^^^ « 1- It is obvious that stronger random potential and narrower band lead 
to a deeper mobility edge. The lower mobility edge of the conduction band can be obtained 
similarly. The energy difference between the lower mobility edge of the conduction band 
and the upper edge of the valence band is 

A™-G- + [i-^ + i^] (5) 

where G*^ is the band gap of LRC. Because the van Hove singularity is smeared out in AS, 
gap in amorphous solid is ambiguous. A^ can be defined in a simulation by identifying two 
edge states. 

In the middle of a band k^a = ^, the group velocity reaches its maximum By Eq. 
(IHI), to localize the states in the middle of the nth band, we need > vr. States in the 

^' — ' ' ' zl„ a ~ 

middle of a band are most difficult to localize. If those states are localized, the whole band 
is localized. A stronger localization condition is Ak ~ k. In the middle of band = f ^, by 
Eq. ([2]) the change in wave vector is It leads to the condition to localize a whole band 



^ > I (smaller than ^ ~ 6 — 34) 13|. The deeper localized states in AS are generated 
by the deeper Bloch states of LRC, are spread in several distorted regions. Because current 
local description only considers the states localized in one distorted region, we cannot expect 
a better estimate. 

The IPR Ij of a localized eigenstate ipj could be approximated asjll] ~ is the 

localization length of ipj. If a Bloch wave ip"^^ suffers a phase shift vr by some distorted 
region to produce ipj, it is localized in range : ^jAk ~ n. The change in wave vector is 



Ak~ 



71 TxVkEk Tczlna sin ka , , 

= -Pl FL 



(~ is obtained under TBA). According to Eq.([T]), F ~ e, e is the relative change in lattice 
constant. To minimize the free energy, a denser region with shorter bonds and small angles 
will gradually decay away toward the mean density rather than exhibit an abrupt transit 
to a diluter region and vice versa. Therefore the size L of a denser distorted region is 
proportional to e. Eq. ([6]) indicates ^ ~ ^j^. The advantage of Eq.Q is that it reveals the 
role of the coordination number z and the transition integral /. The dependence on k (wave 
length and propagation direction of Bloch wave) is also displayed in Eq. ([6]): close to band 
edge of LRC, /ca ~ or ^, the localization length is small and IPR is high (see Fig. [1]). 
Making use of Eqs. and (jl]), 

UEy) = -yj^ll - ( j^^ ) ] (7) 

fej^g is the location of the mobility edge of valence band. When we approach from the 
upper side with higher energy, it is easy to find C,j — * L from Eq. ([7j), localization length 
^ approach to the size L of whole sample as {E^j — &me)"; where | < a < 1, it is close 
the lower bound of previous works jl^. The trend expressed by is consistent with a 



simulation based upon time-dependent Schrodinger equation [15]. 

For a localized state derived from Bloch wave ip^ in LRC, the energy dependence of IPR 
can be found 

;(^,) . i^^^f)' (8) 

This is a new prediction of our work. Eqs. (I7l) and (IE]) are not quite satisfied because Ey 
is the corresponding energy level in LRC, not the eigenvalue of the localized state ipj. It 
can be cured by taking into account energy level shift caused by the disorder in AS relative 

n 

to LRC. Fig. [1] shows IPR vs. eigenvalues in a 512-atom model of a-Sif?^. As expected 
from Eq. ([8]), IPR decreases from highest values at band edges to small values in the band 
interior. The functional form ([8]) fits the simulation rather well. 

According to Eq. (IHl), the least squares fitting parameters in Fig. [T]are {FL)v = 1.256eV, 
{zl)v = 3.185eV, E^ = -7.390eV, {FL)c = 1.437eV, {zl)c = 3.502eV, E^ = -l.OSOeV. 
The width of valence band of c-Si is about 2.7eV, the width of conduction band is about 



2.3eV 



161 ]. The fit parameters are reasonable-something like what we expect for Si. Gap for 




FIG. 1: IPR of 512-atom model of a-Si, dots from ab initio calculation [9[], dashed line and solid 
line are from two parameter (FL and zl) least squares fit and eye guide fit with Eq.Q. 



c-Si is 1.12eV 



16|, using above parameters with help of Eq.dS]), the distance between mobility 



edges is 2.205eV. Result from LRC model falls in the range 1.58-2.43 eV of the observed 
optical gap 17, Q, 3. 



In a distorted region of a-Si where bonds are shortened, valence states have more ampli- 
tude in the middle of bonds. Random potential Va—Vc (the difference between the amorphous 
and crystalline potentials) is important only in the middle of bonds rather than close to the 
core of atoms. Electrons will feel Va — Vc more than a region where bonds more close normal. 
Valence tail states are easily localized in a distorted region with shorter bonds 

On 

the other hand, in a distorted region with longer bonds, the conduction levels are lowered 
and the probability of conduction electrons staying in the middle of nearest neighbor atoms 
becomes larger than a region where bonds are closer to the mean. Conduction tail states 
are more readily localized in a distorted region with longer bonds and large angles [?, l2oj |. 

The effect of three- and four- point correlation on the shape of band tail is subtle: localized 
states adhere to ID filaments in AS network[8]. In the spirit of scattering theory of line 
shape jlll, the decay rate E^^^"^ of valence (conduction) tail can be derived from the relative 
shift of energy levels of LRC. Suppose Ah is the distribution width of bond length (BL), the 
blurring 5k in wave vector k is ^k. The shift of level (E^) for a Bloch state ipl ii'k) 
in valence (conduction) band by the disorder in AS is AEl^'^'' = J dT{Va — V^)|V^fc*'^''P- The 
relative level shift due to this BL distribution is 6k-^AE^. It is easy to see Va — VcCc ^Vc. 
Then 

If we make a correspondence between structural disorder ^IKI and on-site spread W of 



evels, Eq.Q is comparable to Eu ~ 0.5^ {B is the band width) 



221 and Er, ~ ^— ^ 



231], where L and W are correlation length and variance of random potential. Eq. ([9]) is 



also consistent with an assumption of Cody et. al. to explain their absorption edge data 
in a-Si:H|24]. Since the width of BL distribution is ~ 0.1 and \Vc\ ~ 1 — lOeV, the 
order of magnitude of mobility edge should be several tenth eV to leV, so that 

the decay rate of band tails is around several tens to several hundred meV. Both agree with 
experimental observations [2^. Eq.Q indicates E^^'"^ is proportional to static disorder that 
is characterized by (^)^, in consistent with the fact that E^^^^ of a-Si:H increased with 
deposition power |25|]. A6 and h could also be explained as the width and the average value 
of BA distribution. 

Because local compression is compensated by adjacent local tensile in AS, E)^ ~ ^E^, 
where is an order one dimensionless constant characterizing the peak (node) of 

valence (conduction) states. In a-Si and a-Si:H, random potential {Va — K) has larger 
distortion in the middle of Si-Si bonds, since valence states are more in the middle of bonds 
than conduction states 6|, they feel the distortion more. Therefore > q,^ . One expect E^^ 
> E^. This agrees with measurements in a-Si:H: E^ ~43-103meV vs. E^ ~27-37 meV, 
linear relation among Eij and E^j has also been observed [251]. 




FIG. 2: and vs. al^^^: 6 squares are extracted from TBA, dotted line and solid line are 
least square fits with and without (0,0) points. 



To test correctness of 



with 20,000 atoms 



3, 



j ^q.f[9l) , we undertook a TBA calculation for DOS of six a-Si models 
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, l26|. E^^'^\ the width Ccose of BA distributions and the width Ab 
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of BL distribution are extracted. Figl2] clearly shows good linear relation between (EjJ) 
and cr^osg, curves pass origin {E^^'~^^ is zero for crystal) as displayed in Eq. iQ. It can be 
further tested in ion implanted samples, where a continuous increase disorder from crystal 



to amorphous are realized by increasing the dose 



27[. The E}^ (E^) vs. {Abf curve does not 



pass origin (not showing here), this is an indication that BA disorder is a little more decisive 
in determine the shape of a band tail than BL disorder for a well relaxed structure [?. I20I. 



The electron-phonon interaction is strong in AS[28|. At finite temperature, the displace- 
ment of an atom in AS deviate from the position in the LRC at zero temperature is a 
vector sum of the static displacement u<j and thermal vibration displacement UT{t) from 
the zero temperature configuration of AS, t is the time moment. In ordinary absorption 
experiment, time interval T is much longer than the period of the slowest mode, therefore 
Efj = ^ (jt<;''"( "''"^"^'-*^ )^|V^|. Atoms vibrate around their equilibrium points in AS, the 
time average of the cross term ■ U2-(t) is zero. Thus Urbach energy from static disorder 
and from thermal disorder is additive[24] E^ = E^^ + E^j,, E^^ = <;'^(^)^| K]. Thermal 
part E^r,. = <^'"^|Vc|, = ^ dt[uT{t)]'^ is the long time average of the square of ampli- 
tude of vibration. An ultra-fast probe of absorption edge may find oscillating in E^. Since 
^ a ^a^H, Be is binding energy in the diluter regions where conduction tail states are 
localized, E^rp = q'-^ksT^-^. E^ linearly increases with temperature. Similarly result holds 
for EK. The is consistent with the fact that above 350K absorption edge linearly increase 



with ksT in a-Si: H[29|,l30fl. Because By > Be, E\j is more susceptible to thermal disorder 
than Eij 3l|, as observed in ref. 30 1. 

For realistic amorphous solid with topological disorder, by viewing an AS as many dis- 
torted regions relative to corresponding LRC, we push forward essential understanding on 
localized states confined in one distorted region. The predicted IPR, mobility edge, the 
dependence on static disorder and on temperature of the decay rate of band tails agree with 
available experiments and simulations. We explained the fact that valence tail states are 
more localized in a denser region with smaller BA and shorter BL and conduction tail states 
a.e ,„o..e loca..ed i„ d.K.te. .eg,o. witK .o^g. BL a.d l.,e. BA . a-SiQ Q. Local- 
ized states in several distorted regions and other problems involving global topology will be 
addressed in future. 
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